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ABSTRACT 

We report on phase-referenced VLBI radio observations of the Type lib supernova 2011dh, at times 
t = 83 days and 179 days after the explosion and at frequencies, respectively, of 22.2 and 8.4 GHz. We 
detected SN 2011dh at both epochs. At the first epoch only an upper limit on SN 2011dh's angular size 
was obtained, but at the second epoch, we determine the angular radius SN 2011dh's radio emission 
to be 0.25 ± 0.08 mas by fitting a spherical shell model directly to the visibility measurements. At a 
distance of 8.4 Mpc this angular radius corresponds to a time-averaged (since t = 0) expansion velocity 
of the forward shock of 21000 ±7000 km s~F Our measured values of the radius of the emission region 
are in excellent agreement with those derived from fitting synchrotron self-absorbed models to the 
radio spectral energy distribution, providing strong confirmation for the latter method of estimating 
the radius. We find that SN 2011dh's radius evolves in a power-law fashion, with R cx i°- 92±01 °. 

Subject headings: supernovae: individual (SN2011dh) — radio continuum: general 



1. INTRODUCTION 

Supernova SN 2011dh was discovered on 2011 May 31 
by the amat eur astronomer A madee Riou in the nearby 
galaxy M51 (iGriga et al.ll2011[). whic h is at a distance of 
8.4 ± 0.6 Mpc (jFeldmeier et alJll997t IVinko et al.ll2012h . 
The supernova was soon confirmed u sing pre- and post- 
discovery imaging (Gr iga et al.ll20lil ). In particular the 
optical transie nt was also detec ted by the Palomar Tran- 
sient Factory (|Law et al.1 120091 ) shortly after the initial 
discovery. The explosion date is tightly cons trained to be 
betw een 2011 May 31.275 and 31.893 UT (jArcavi et al.1 
2011). We will adopt the (rounded) midpoint of this 
interval of to = May 31.6 UT. The supernova was coin- 
cident with the eastern spiral arm of M51, and had an 
estimated apparent magnitude of ~ 14 mag (unfiltered) 
implying an absolute magnitude of roughly ~ — 16 mag. 
We note that M51 is also host to an earlier supernova, 
SN 19941, which however was of Type Ib/c. 

Initially, S N 201 ldh was spectrq scopically classified 
as Type IIP (|Silverman et al.|[2011l) . but further spec- 
troscopy, which showed helium absorpt ion features, 
caused a re-classification as a Type lib (jArcavi et al.1 
l20Tfl : M anon eTaDHoiil). A maximum expansion ve- 
locity of ^17000 k m s -1 was estimated from the blue 
edge of the Ha line ([Silverman et al.lfeoill : lArcavi et al.1 

1ml). 

Radio emission was detected on June 4, only a few 
days after the exp losion, at centimeter wavelengths 
(|Horesh et al.l l2011af) with the National Radio Astron- 
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omy Obs ervator J] (NRAO) Expanded Very Large Array 
EVLA; | Perlev et alj 1201 If) , at millimeter wavelengths 
Horesh et al.ll2011bf ) using the Combined Array for Re- 
search in Millimeter- w ave Astronomy (CARM A), and at 
sub-mm wavelengths ([Soderberg et al.l 120111 ) using the 
Submillimeter Array (SMA). The initia l radio and mm- 
band observations were presented in ISoderberg et al.1 
(|2011f) . Further radio flux-density measurements as well 
as modeling of the light curve a re presented i n a co m- 
panion paper to the present one. iKrauss et all ((2012) . 

Although a yellow supergiant star was identified on 
pre-explosion Hubble Space Telescope (HST) images as 
a pos sible progenitor (IVan Dvk et al.|[20TTt IMaund et al.1 
1201 If ) , subsequent wor k bvlArcavi et alj (l201lh suggested 
a smaller progenitor. ISoderberg et a l. (2011) examined 
the early data at radio through X-ray wavelengths and 
concluded that these also suggest a small progenitor star 
with a stellar radius of ~ 10 11 cm and were not consistent 
with a radius of order 10 13 cm as would be expected for 
a yellow supergiant. The yellow supergiant identified in 
the HST images would then have to be either a binary 
companion of, or possibly a chance superposition with, 
t he true progenitor star 

IChevalier fc Soderberg! (|2010f ) proposed that Type lib 
supernovae (SNe lib) are divided into two distinct sub- 
classes. The first, SN cllb, have compact progeni- 
tors, with stellar radii ~ 1 11 cm, (e.g., SN 20 01ig, 
iRvder et alj[200l SN 2003bg.lSoderberg et alj2006t and 
SN 2008ax. lRoming et al.ll2009D . Thev are characterized 
by high shock velocities of ~ 0.1c, and radio light curves 
which show deviations from a power-law decay at late 
times . The seco nd, SN ellb, (e.g., SN 1993 J,|Ba rtel et al.1 
[20021 SN 2001gd |Perez-Torres et aT]2005l) have extended 
progenitors, with stellar radii ~ 10 13 cm, and are charac- 
terized by somewhat slower shock velocities and smooth 
radio light curves. 

The size and expansion velocity of the shockfront is a 

7 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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basic characteristic distinguishing different supernovae, 
and it is therefore important to determine it observa- 
tionally as directly as possible. In particular, to deter- 
mine whether SN 2011dh conforms to the characteristics 
of Chevalier & Soderberg's proposed Type SN cllb class 
requires measuring the expansion velocity of the shock 
front. Very long-baseline- interferometry (VLBI) obser- 
vations are the most direct way of making this measure- 
ment (see e.g., IBietenholz et al.l 120101 : iBrunthaler et al.l 
120101: IBietenholz! 120081 ). Unlike the optical emission, 
which mostly originates in the denser and more slowly 
moving inner ejecta, the radio emission generally traces 
the fastest ejecta. The radio emission is thought to 
originate in the region between the forward and reverse 
shocks. In the part icularly well-studied case of SN 1993 J, 
iBartel et al.l ()2007f ) show that there is a close relationship 
between the outer boundary of the radio emission and the 
location of the forward shock. 

However, even with the high resolution afforded by 
VLBI, it is generally difficult to measure the radius 
of the radio-emitting region, and thus of the forward 
shock, especially early on in the evolution of a super- 
nova when its size is still sm all. Chevalier (1998; see also 
iChevalier fc Fra nsson 2006) has shown that, for a super- 
nova spectrum dominated by synchrotron self- absorption 
(SSA), the radius of the emitting region, i?, (assumed 
spherical), as well as the magnetic field can be deter- 
mined from two observable quantities: the frequency and 
spectral luminosity of the peak in the radio spectral en- 
ergy distribution (SED). The observables require only 
measurements of the total flux density at different fre- 
quencies, and thus are not dependent on spatially resolv- 
ing the supernova. In the case that significant free-free 
absorption is present in addition to SSA, a lower limit 
on the shockfront radius is obtained. 

This calculation of the radius is fairly robust. However, 
in addition to the mentioned observables and the dis- 
tance, D, which is usually fairly well constrained, the cal- 
culation does also rely on assumed values for two poorly 
known parameters. The first of these, /, is the filling 
factor of the radio emission, while the seconcQ, tp, is the 
ratio of relativistic electron energy density to magnetic 
energy density in the post-shoc k region. 

In the case of SN 2011dh, ISoderberg et~aTl (l20ll 
showed that the cm and mm-wave radio flux densities 
are consistent with a synchrotron self-absorbed spec- 
trum at t ~ 5 and 17 days, and determined that R ~ 
3.7 x 10 15 cm at the latt er epoch. In our companion pa- 
per, iKrauss et al.l (|2012l ) continue this work through to 
i = 92 days, and present detailed EVLA monitoring of 
the multi-frequency radio light curves, as well as mod- 
eling the light curves to show that the evolution of R. 
as calculated from the radio SED, is consistent with a 
power-law evolution with R oc t 9 . 

We report in this paper on VLBI observations of 
SN 201 ldh obtained at t = 83 and 179 days after the 
explosion, which have allowed us to directly measure (or 
at least constrain) the angular size of SN 201 ldh at those 
epochs. Early results from t he first set of obs e rvation s 
have already been reported in IBietenholz et al.l (|2011b| ) . 

8 Note that IKrauss et alj <20l3) . ISoderberg et al.l lr201ll) and 
ICheval ier ( 1998) use the symbol a for this parameter, which symbol 
we use here to denote the radio spectral index. 



2. RADIO LIGHT CURVES 

To set the context for, and aid in the interpretation of, 
the VLBI results, we present in Figure [TJ the radio light 
curves of SN 2011dh at 22.2 and 8.4 GHz, the two fre- 
quencies at which we carried out the VLBI observations. 
The flux density values at 22.2 GHz were logarithmically 
interpolated in frequenc y between 20 . 5 and 25.0 GHz 
from those presented in IKrauss et al.l (|2012f l. with the 
exception of the last one at t = 183 days. For that last 
one, we obtained additional EVLA observations on 2011 
November 30 at 20.5 and 25 GHz. The data redu ction 
procedure was as described in IKrauss et al.l (|2012t ). We 
obtained flux densities of of 722 ± 72 /zJy at 20.5 GHz 
and 551 ± 55 /iJy at 25 GHz, with the uncertainties be- 
ing intended as standard errors and dominated by an 
estimated 10% systematic contribution. The value for 
the flux density at the time of the 22.2-GHz VLBI ob- 
servations at t = 92 days was then interpolated logarith- 
mically first in frequency and then in time between the 
adjacent values. 

The 8.4-GHz m easurements were also taken from 
IKrauss et~ al. (2012). Near the time of the second VLBI 
observations, the VLA was in the compact D array con- 
figuration, with resolution insufficient to overcome con- 
fusion at 8.4 GHz, so a direct measurement could not be 
obtained. We therefore estimated the 8.4-GHz flux den- 
sity by first interpolating the 2 2. 2- GHz values logarith- 
mically in time to t = 179 days, and then extrapolating 
to 8.4 GHz by assuming the same radio spectral index of 
a = — 1 (where S v oc u a ) as was found appropriate for 
the opti cally thin part of t he spectrum through to t = 
92 days (jKrauss et all2012T ). Although the optically-thin 
spectral index h as been seen to vary with time in some 
supernovae (e.g. . IBartel et al.ll2002t IBartel fc Bietenholzl 
2008), large variations are not expected, so the re is no 
reason to think the model of lKrauss et al.l (|2012f) can not 
be extrapolated to t — 179 days, so our estimate of the 
8.4-GHz flux density at that time is of sufficient accu- 
racy for comparing to the flux density recovered from 
the VLBI observations. 

3. VLBI OBSERVATIONS 

Our two sessions of VLBI imaging observations of 
SN 201 ldh were carried out using the high-sensitivity 
array, which consisted of the NRAO Very Long Baseline 
Array (VLBA; 10 x 25-m diameter, distributed across 
the United States), the NRAO Robert C. Byrd (-105 m 
diameter) telescope and the Effelsberg (100 m diameter) 
telescope. Both sessions were 13 hr in length, with the 
first using an observing frequency of 22.2 GHz on 2011 
August 21 UTC, and the second using an observing fre- 
quency of 8.4 GHz on 2011 November 26 UTC. The age 
of the supernova was 83 and 179 days, respectively, at 
the midpoints of our two sessions. 

In each set of observations, we included several "geode- 
tic blocks," where we observed 12-15 bright sources from 
the Interna tional Celestial Reference Frame (ICRF) list 
of sources (|Fev et al. 2009) over a period of —45 min- 
utes to measure the troposph eric zenith delay and 
clock offsets at each ante nna (see Brunth aler et alj|2005t 
iReid fc Brunthalerll2004f) . In these geodetic blocks we 
recorded 8 intermediate frequencies (IFs) with frequency 
centers spread over a ~400 MHz range, with each IF 
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Figure 1. Radio light curves of SN 2011dh at 22.2 and 8.4 GHz, as obtained from EVLA observations. Flux density measurements are 
shown by solid circles and squares. See text for details of how t he flu x density values were obtained. With the exception of the one at age 
= 183 days, they were taken or interpolated from Krauss ct al. (2012). The open square represents the interpolated value at our first VLBI 
observing date and frequency (22.2 GHz at age = 83 days) while the open circle and dashed line show the extrapolation to our second 
VLBI session (8.4 GHz at age = 179 days). The age of the supernova is calculated from the explosion date of 2011 May 31.6 UT. 



covering a bandwidth of 8 MHz. 

For the observations of SN 201 ldh we recorded a con- 
tiguous bandwidth of 64 MHz in each of the two senses of 
circular polarization with two-bit sampling, for a total 
bit rate of 512 Mbit s" 1 . We used J1332+4722 (ICRF 
J133245. 2+472222), 0?5 away from SN 2011dh, as a pri- 
mary phase cal ibrator. This source is an ICRF source 
(|Fev et al.l 2009) with a position known to —70 /ias. Any 
positions in this paper are calculated by taking the posi- 
tion of J1332+4722 to be RA = 13 h 32 m 45?24642, deck 
= 47° 22' 22"6670 for J1332+4722. For the first session, 
at 22 GHz, we used a cycle time of —110 s, with —60 s du- 
ration on SN 2011dh and -50 s on J1332+4722, while for 
the second session at 8.4 GHz, we used a somewhat longer 
cycle time of -170 s, with -120 s spent on SN 2011dh. 

We considered as a possible phas e-reference source th e 
weak nuclear radio source of M51 (IMaddox et al J [20071. 
which is only —3' from SN 201 ldh. This nuclear source, 
however, has a low 5- GHz peak brightness of <1 mJy, 
as well as having a steep spectral index, and is therefore 
too weak to use for phase referencing at our observing 
frequencies of 8.4 and 22.2 GHz. 

During each of the observing runs, we also spent three 
— 20-m periods observing an astrometric "check" source, 

9 We note that parallel recordings were made using the higher 
bandwidth of 256 MHz per polarization, corresponding to a bit 
rate of 2 Gbit s — 1 , with the experimental new Mark 5C recording 
units. Various technical difficulties with these recordings, however, 
rendered the correlated data unreliable, and we base the results in 
this paper only on the data recorded at 512 Mbit s — 1 using the 
well-tried Mark 5B recording systems. 



the quasar JVAS J 1324+4743, which was about 1?5 away 
from our primary phase-reference source J1332+4722. 
The purpose of the observations of J 1324+4743 was to 
check the quality of the phase-referencing and also to pro- 
vide a second astrometric reference source. This check 
source was observed using a similar phase-referencing 
pattern as we used for SN 201 ldh. 

T he VLBI data were correlated with the DiFX correla- 
tor (jDeller et al.ll20lil ) , and the analysis carried out with 
NRAO's Astronomi cal Image Processing System (AIPS) 
and ParselToungue ijKettenis et a l. 2006). We calibrated 
both sets of observations using standard procedures, 
making a correction for the dispersive ionospheric delay 
using the AIPS task TECOR, and solving for the zenith 
tropospheric delay on the basis of our geodetic-block ob- 
servations. We discarded any SN 2011dh visibility data 
obtained when either of two the antennas involved was 
observing at elevations below 10°. 

The initial flux density calibration was done through 
measurements of the system temperature at each tele- 
scope, and then improved through self-calibration of the 
primary reference source J1332+4722. This source is 
slightly resolved, as can be seen on the images in the 
VLBA calibrator list data-bas43, where a weak exten- 
sion or second component is visible —2 mas to the west- 
southwest of the peak. We see a similar structure in the 
images made from our data at both 22 and 8.4 GHz. Our 
final amplitude and phase calibration at both frequen- 
cies was derived using a CLEAN model of this source, 

10 http://www.vlba.nrao.edu/astro/calib 
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with the peak-brightness point in the image being placed 
at the nominal coordinates given abovePl. Finally this 
calibration was interpolated to the intervening scans of 
SN 2011dh. 

4. RESULTS 

4.1. Results from the 2011 August 24 Observations at 

22 GHz 

Unfortunately, the weather was generally poor for our 
first set of observations, at 22 GHz, resulting in higher 
than usual noise levels: the rms background bright- 
ness in a naturally- weighted image was 92 /iJy beam -1 . 
Nonetheless SN 2011dh was clearly detected, with a peak 
brightness of 630 /iJy beam -1 , which was 6.8 times the 
rms background. As mentio ned above, early results f rom 
this epoch were reported in IBietenholz et al.l (|2011b[ ) . 

For marginally resolved sources, such as SN 201 ldh, 
the best values for the source size and VLBI flux density 
come from fitting models directly to the visibility data, 
rather than from imaging. We chose as a model the pro- 
jection of an optically-thin spherical shell of uniform vol- 
ume emissivity, with an outer radius of 1.25 x the inner 
radiutQ Such a model has been foun d to be appropri- 
ate for other radio supernovae ( see e.g..[Bictcnhol z~et al.l 
12003 iBartel fc Bietenholzll2008l ). The Fourier transform 
of this shell model is then fitted to the visibility measure- 
ments by least squares. For a partially resolved source 
such as SN 201 ldh, the exact model geometry is not criti- 
cal, and our shell model will give a reasonable estimate of 
the size of any circularly symmetric source, with a scaling 
factor of order un ity dependent on t he exact morphology 
(see discussion in IBartel et al.ll2~002l ). 

For the observations of SN 2011dh at t = 83 days, 
our best-fit spherical shell model had a total flux density 
of 650 ± 140 /iJy (statistical and systematic uncertainty 
combined). We note that the total flux density as inter- 
polated from the VLA measurements was 1100 /xJy (see 
Section[2]and Figure[T]), with an estimated uncertainty of 
10%, and is therefore higher than that recovered from the 
VLBI observations by a combined 2.5er. This discrepancy 
suggests that there is likely some loss of phase-coherence 
in the VLBI observations, which is not unexpected given 
the relatively poor weather. 

We obtained a value of O.ll^g'Jg mas for the outer 
angular radius of SN 201 ldh at this epoch, where the 
listed uncertainty is intended as a standard error, and 
consists of both the statistical and systematic compo- 
nents added in quadrature. We estimated the systematic 
component as follows: the largest systematic error most 
likely arises from the uncertainty in the antenna ampli- 
tude gains, which are only imprecisely known but are 
correlated with the fitted source size for marginally re- 
solved sources, as well as from the likely coherence losses 

11 The deviation of the source geometry of J1332+4722 from 
a point source is small enough so that the effect of using a point 
model in the solutions for delay and delay rate made using FRING 
solutions is negligible. 

12 Our results do not depend significantly on the assumed ratio 
between inner and outer radii, as the effect of reasonable varia- 
tions in this ratio is considerably less than our stated uncertainties. 
For a discussion of uncertainties on the shell-size obtained through 
u-v plane modclfitting compared with those obtained in the image 
plane for the case of SN 1993 J, showing that superior results are 
obtained using the former, sec Bietenholz et al. (2011a). 



mentioned above. Although the effect of the latter is 
difficult to estimate reliably, the first order effect would 
be a reduction in the visibility amplitude. We accord- 
ingly estimated the effects of both amplitude-gain errors 
and coherence loss from a small Monte- Carlo simulation 
where we randomly varied the amplitude gains of the 
individual telescopes by 20% rms, to arrive at the uncer- 
tainty in the angular radius listed above. The supernova 
is therefore not significantly resolved at this epoch, and 
the 3cr upper limit on the outer radius was 0.38 mas. 

We also take the fitted center position of the model 
as our best estimate of the center-position of SN 201 ldh. 
For this epoch, this position was RA = 13 h 30 m 5! 105548 
and deck = 47° 10' 10"92273. The statistical uncertainty 
on this position is ~50 /i&s, however, a realistic position 
uncertainty will have substantial contributions from un- 
modeled effects of the troposphere, ionosphere, station 
position errors, and possible evolution of the reference 
sources. We estimate the total uncertainty in the po- 
sition relative to that of J1332+4722 to be ^70 //as. 
Our measured position therefore agrees to within less 
than the combined unce rtainties with that measured by 
iMarti'-Vidal et al.l ()2011h at t = 14 days. 

4.2. Results from the 2011 November 26 Observations 
at 8.4 GHz 

In the second set of VLBI observations at 8.4 GHz, 
SN 201 ldh was detected with good signal-to- noise ra- 
tio. We show the VLBI image of SN 201 ldh in Fig- 
ure H The fitted FWHM size of the CLEAN beam 
was 0.72 x 0.51 mas at p. a. —30°. For presentation 
we slightly super-resolved the image by restoring with 
a round clean beam of FWHM 0.51 mas. We em- 
phasize that our measurement of the angular size of 
SN 201 ldh, as presented below, is entirely independent 
of this super-resolution. The peak brightness in this im- 
age was 1160 /iJy beam -1 , while the rms background 
brightness was 22 /iJy beam -1 . 

Since SN 201 ldh is only partly resolved, we again turn 
to fitting models directly to the visibility measurements 
in order to accurately measure the size of the source. Due 
to the good signal-to-noise ratio obtained at this epoch, 
the source size could be fairly accurately constrained by 
this method. In the model-fitting, we used the square- 
root of the uncertainties of the visibility measurements, 
which procedure, though statistically less efficient, in- 
creases the robustness of results in the likely case that 
the errors in the visibility measurements are not purely 
statistical. We again used a spherical shell model with 
an outer radius 1.25 times the inner radius. 

From the model-fitting, we obtained a total flux den- 
sity of 1540 ± 190 /iJy (an almost identical value was ob- 
tained from imaging) . This value is only 7% below that of 
1.64 mjy extrapolated for this epoch and frequency from 
the EVLA measurements (see Section [5] and Figure [T]), 
although we note that the extrapolation EVLA value is 
somewhat uncertain in that it depends on the errors in 
the individual flux density measurements as well as the 
assumption that the radio spectral index does not vary 
in time. We find therefore that the flux density recovered 
from the VLBI measurements is in reasonable agreement 
with that extrapolated from the EVLA measurements, 
and that there is no reason to suspect significant corre- 
lation losses. 
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Figure 2. The CLEAN image of SN 2011dh made from VLBI observations at 8.4 GHz on 2011 November 26 UTC. To obtain more robust 
imaging, we started with the square root of the visibility weights and then used uniform weighting. The fitted CLEAN beam FWHM was 
0.72 X 0.51 mas at p. a. —30°, and we super-resolved slightly by using a round CLEAN beam with FWHM of 0.51 mas, indicated at lower 
left. The contours are drawn at 10, 20, 30, 50, 70 and 90% of the peak brightness, with the 50% contour being emphasized. The peak 
brightness was 1160 fijy beam -1 , and the background rms brightness was 22 /jjy beam -1 . The greyscale is labeled in ^tjy beam -1 . North 
is up and east is to the left, and the coordinate origin is at the fitted center position of SN 2011dh, which was at RA 13 h 30 m 05?105548, 
decl. 47° 10' 10'.'92273. 

referenced data, we found that the fitted FWHM major 
axis was 25% larger than when we fit to the phase self- 
calibrated data. This implies that, at least in the case 
of our check source J1324+4743, the coherence losses 
due to phase-referencing can result in an error of 25% 
on the source size. We note however, that J1324+4743 
was almost three times farther from the phase-reference 
source (J1332+4722) than was SN 201 ldh. Further- 
more, the u-v coverage for J1324+4743 was much poorer 
since we only observed it for three brief periods. We 
expect therefore that the additional errors in the size of 
SN 201 ldh due to coherence loss are considerably smaller 
than 25%, and thus considerably smaller than the un- 
certainties from other causes discussed above. We do 
not, therefore, find any reason to suspect that coherence 
losses due to poor phase-referencing would substantially 
increase the uncertainty we give above for the angular 
size of SN 201 ldh. 

5. DISCUSSION 

We obtained VLBI observations of the nearby Type lib 
supernova 201 ldh, with the primary goal of obtaining a 
direct observational constraint on the expansion speed by 
measuring the angular size of SN 2011dh. VLBI observa- 
tions are crucial, as they are the only means to directly 



We again also obtained a center position of SN 2011dh 
from this fit, which was consistent within the uncertain- 
ties with that obtained for our first ep och as well as with 
the po sition obtained at t — 14 days bv lMarti'-Vidal et"al"1 
(|2011| ). The formal difference in position over the 165- 
day interval between the Marti- Vidal ct al. observation 
and the present one was 110 ± 125 /ias. 

Finally, we obtained an outer radius for SN 201 ldh of 
0.25 ± 0.08 mas, with the uncertainties again including 
both statistical and systematic contributions and derived 
as in the previous section. At a distance of 8.4 Mpc, 
this angular radius corresponds to an average expansion 
velocity since the explosion of 21000 ± 7000 km s _1 . 

To determine the quality of the phase referencing, 
upon which the above determinations of the position 
and angular radius depend, we examined the obser- 
vations of our check source, J1324+4743. From the 
strictly phase-referenced data for the check source we 
found a total flux density of 101 mJy, while for the 
phase self-calibrated data we found a 36% higher value 
of 137 mJy. We fitted an elliptical Gaussian model, ap- 
propriate for a marginally resolved QSO, directly to the 
visibilities in a similar fashion as the shell model was fit 
to the SN 2011dh data. When fitting the strictly phase- 
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measure the size, expansion speed and perhaps the ge- 
ometry of the radio emission region. We obtained an 
upper limit on the radius of SN 201 ldh at t = 83 days 
and a measurement of the radius at t = 179 days. We 
plot these two values as a function of time in Figure |3] 
(red squares). 

As mentioned in the introduction, the radius of a su- 
pernova can also be calculated from broadband, but not 
spatially-resolved, observations of the radio SED (spec- 
tral energy distribution) , under the assumption of a spec- 
trum dominated by SSA (synchrotron self-absorption). 
In our companion pa per Krauss et al. (2012; see also 
ISoderberg etalll20ll we used such observations of the 
SED to calculate the shockfront radius for SN 2011dh 
through to t = 92. In particular, the radius is calculated 
follow ing Chevalier & Fransson (2006; see also iChevalierl 
1998): 




where S v is the observed flux density at the peak of the 
synchrotron spectrum, which occurs at frequency v op , D 
is the distance, and / and ip are, as mentioned above, 
the filling factor and the ratio of the energy density in 
relativistic particles to that of the magnetic field, respec- 
tively. 

For / we have some observational constraints at least in 
the well-studied case of another SN lib, SN 1993J, where 
the well-resolved VLBI images suggest a very spherical 
shell structure with a shell thickness of 20% to 33% of the 
outer radius (e.g.. IBartel et all 120021 : IBietenholz et al.l 
I2011at iMarti'-Vidal et alj|2010f) Such a geometry would 
have a gross value of / in the range of 49% to 70%, al- 
though it is possible that small scale dumpiness within 
this overall geometry could reduce / below those values. 
For ip, the commonly used value is ip — 1, correspond- 
ing to equipartition . Here, and in our companion paper 
(|Krauss et alj |2012) we therefore take / = 0.5 and if: = 1. 

The radius of SN 2011dh's forward shock was calcu- 
lated from the SED as obtained from EVLA and SMA 
observations using this m ethod at various times be- 
tween t = 4 and 92 days by ISoderberg et al.l (|2011D and 
iKrauss et al.l (|2012l ). We plot these values also in Fig- 
ure |31 (green circles). 

Assuming that the supernova ejecta and circumstellar 
material have power-law density profiles, then evolution 
of the forward shock radius with time follows a power- law 
form and can be expressed simply as R cx i m , where m is 
known as the deceleration parameter. If the circumstellar 
density is, as expected for a stellar wind, oc R~ 2 , and the 
outer portions of the ejecta are also characterized by a 
power-law density distribution with in dex n, then it ca n 
be shown that m = (n — 3)/(rt - 2) (|Chevaliedll982T) . 
Indeed, we found in IKrauss et al.l (|2012l ). that at least 
up to t = 92 days the evolution of SN 2011dh's radius 
was consistent with m = 0.9, implying that n a 12. 

In order to determine m from our various measure- 
ments of SN 2011dh's radius, we performed a weighted 
least-squares fit to the complete set of values for i?, in- 



cluding both those determined from our VLBI measure- 
ments and those from fitting the SED. The fitted line 
is also plotted in Figure [3] (in blue) . As can be seen in 
Figure El the expansion of SN 201 ldh is in fact well de- 
scribed by a power-law. The fit gives to = 0.92 ± 0.03 
and a radius at t = 30 days of (5.85 ± 0.13) x 10 15 cm, 
where the uncertainties are statistical ones derived from 
the scatter of the fit. 

In addition to the power-law nature of the expansion, 
it can be seen from Figure [3] that there is excellent agree- 
ment between the radii measured with VLBI with those 
calculated from the radio SED under the assumption of 
SSA. This agreement provides strong validation for the 
latter method of calculating the shockfront radius. In- 
deed, SN 2011dh represents so far the best example for 
directly comparing the radii of the shock wave deter- 
mined in these two different fashions. 

The radii derived from fits to the SED are uncertain 
due to a number of systematic factors. Firstly, they 
are calculated from the fitting of a model SSA spec- 
trum to the observed flux density values. In the case 
of SN 201 ldh, although the SSA model is an excellent 
fit o yer a wide spectral range, the fit is not exact (see 
e.g.. IKrauss et a 1.1 12012;. resulting in estimated system- 
atic uncertainties of 15% and 5%, respectively, on the 
values of v op and S Vop , which result in a ~15% system- 
atic uncertainty on the calculated radius. In Figure|3l the 
plotted error bars, intended as standard errors, include 
this systematic component on the values of R calculated 
from the SED. We performed a Monte-Carlo simulation 
where we varied all the radius measurements according 
to their standard errors, and found an rms scatter of the 
derived values of to of 0.10. It can be seen, however, that 
any such systematic errors do not seem time-dependent, 
in other words likely have only a small effect on the de- 
rived value of to. We therefore conservatively estimate 
the systematic la uncertainty on our fitted value of to 
at 0.10. 

The radii calculated from the SED also depend, as 
noted above, on the poorly known para meter ip, for which 
we too k the equipartition value of unity. ISoderberg et al.l 
(|2011[ ). however, found that for SN 201 ldh the X-ray 
measurements, in conjunction with those in the radio, 
suggested deviations from equipartition, with ip ~ 30. 
Can we constrain ip by comparing the values of the R 
as calculated from the SED with those from VLBI which 
are independent of "0? 

As can be seen from Equation Q] above, the former val- 
ues of R in fact depend only weakly on ip. In particu- 
lar, the change of ip fro m our assumed va l ue of unity to 
ip = 30, as suggested by ISoderberg et all ()201lD . would 
decrease the calculated values of R by 16%, producing 
only a small change in Figured Using such slightly lower 
values of the R determined from the SED (but with the 
original values of R from VLBI) does not change the fit- 
ted value of m by more than the statistical uncertainty, 
and results in a slightly, but not significantly, poorer fit 
of the power-law expansion. We therefore conclude that 
the measurements are not yet of sufficient accuracy to 
usefully constrain ip!. 

As we have shown, the deceleration parameter, to, for 
SN 2011dh seems robustly determined at m ~ 0.9, which 
implies only slightly decelerated expansion, and conse- 
quently a relatively steep radial density profile in the 
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Figure 3. The shock front radii of SN 2011dh as estimated by two independent methods. In red, we plot the values derived from fitting 
spherical shell models to the VLBI visibility-data from this paper. The plotted la error bars include statistical and systematic contributions 
(see text Section[3]for details). In green, we plot the values calculated from the radio spec tral energy distribution unde r the assumption tha t 
it is dominated by SSA (synchrotron self-absorption), with the radius values taken from Sodcrbcrg ct al. (2011) and Krauss ct al. (2012), 
again with uncertainties including statistical and systematic components. The blue line represents our power-law fit to all the values, with 
R = 5.85 X 10 15 (t/30days) 92 cm (see text). 



oute r e.jecta, with po wer-law index n ~ 12. We note 
that lChevalie"r1 ()1982[ ) found that for self-similar models 
with m = 0.9 and n — 12 the ratio of outer shock to in- 
ner shock radii was 1.24, almost identical to our adopted 
value of 1.25. 

The deceleration of SN 201 ldh may show a pattern 
similar to that seen for SN 1993 J, w here m 9 was 
found for approximately the first year (jBartel et al.ll200l: 
IMarti-Vidal et all 12010). albeit with somewhat smaller 
velocities than SN 2011dh. We note that the expansion 
velocities SN 2011dh's shock front, as determined from 
the radio emission, are approximately twice as large as 
the highest velocity of ~ 17000 km s _1 observed in the 
opti cal spectrum (from the blue edge of the Ha emission 
line; iSilverman et al.|[201lt lArcavi et all 1201 ID . This is 
not unexpected as the high- velocity shocked H i s likely to 
be non-radiative (jChevalier fc Soderberd[20Tol) . and the 
e.jecta should be characterized by a steep densi ty profile 
(jBerger et al.l[200l[Chevalier fc Franssonll2006f ). 

Our observations also constrain the proper motion of 
the geometrical center of SN 2011dh: be tween the EVN 
observations of IMarti-Vidal et al.l (|2011l) at t = 14 days 
and our observations at t = 179 days the observed proper 
motion corresponds to a velocity of 9600 ± 10500 km s~ x , 
not significantly different from zero. We note that large 
proper motions of the center of the emission region are 
not expected: the best-determined case is SN 1993J, 
where a peculiar motion of only 320 ± 160 km s _1 was 
observed dBietenholz et al.ll200lD. 

As mentioned, IChevalier fc Soderberd (pOlOt ) have sug- 
gested that type lib supernovae are divided into two 



categories according to whether the progenitor was ex- 
tended or c ompact. Since SN 2011dh had a compac t 
progenitor (jSoderberg et al.l 120111: Krauss et ail 12012ft , 
it would belong to th e latte r category. According to 
IChevalier fc Soderbergi (|2010Q this cate gory is character- 
ized by a high expansion velocity as well as by varia- 
tions in the light curve, and presumably also the decel- 
eration, due to episodic mass-loss from the progenitor. 
SN 2011dh's relatively high expansion velocity has al- 
ready been noted, and is confirmed by our VLBI ob- 
servations. It will, however, be important to continue 
observations of SN 2011dh's radio SED as well as VLBI 
imaging for as long as possible to determine whether the 
remainder of the predictions are borne out, and to in- 
crease our sample of supernovae with observationally de- 
termined expansion curves as well as radio and X-ray 
light curves, since supernovae as nearby and as radio- 
bright as SN 201 ldh are rare events. 
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